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The structure of finite-area topological defects in graphene is described in terms of both the direct honeycomb lattice and its 
dual triangular lattice. Such defects are equivalent to cutting out a patch of graphene and replacing it with a different patch 
with the same number of dangling bonds. An important subset of these defects, bound by a closed loop of alternating 5- and 
7-membered carbon rings, explains most finite- area topological defects that have been experimentally observed. Previously 
unidentified defects seen in scanning tunneling microscope (STM) images of graphene grown on SiC are identified as isolated 

di vacancies or divacancy clusters. 



Producing commercial graphene-based devices will re- 
quire the ability to grow large sheets of high-quality 
graphene. Several techniques for producing graphene 
exist, including mechanical exfoliation from graphite pQ, 
chemical exfoliation from graphite [2 , chemical reduc- 
tion of graphene oxide [3 , segregation of carbon from 
metal crystals [4 , chemical vapor deposition of C onto 
metal surfaces [5j [6], and thermal desorption of Si from 
SiC [3 [8]. Graphene produced by the above methods is 
often found to contain defects, such as vacancies [9] or 
grain boundaries [TQHT2] . Defects decrease the high mo- 
bility of graphene [13]; therefore, it is desirable to reduce 
or eliminate the number of defects. Conversely, one may 
desire to tune the properties of graphitic materials, e.g. 
the bandgap, by intentionally creating and manipulating 
defect structures [T4HT7] . Toward either end, it is neces- 
sary to classify the kinds of defects that form in graphene 
and to correlate them with the growth conditions. 

One class of defects, frequently observed in scanning 
tunneling microscope (STM) images of ultrahigh- vacuum 
graphene growth via Si desorption from SiCpQ [18 , ex- 
hibits regions, roughly several nm in size, of strongly per- 
turbed electronic structure completely surrounded by or- 
dered graphene. The finite range of the electronic struc- 
ture perturbation suggests that these defects could be 
created or healed by the motion of a relatively small 
numbers of atoms, and may therefore be among the most 
important defects in graphene. We recently identified [19] 
the sixfold symmetric defect seen in Ref.[7]as the "flower" 
defect [20, 21 (shown below), a topological defect that 
can be described as the rotation of 24 central atoms in 
ideal graphene by 30°. Other finite-area defects, with 
twofold, and threefold symmetry, are seen in Ref. [7] but 
their structures have not been previously identified [18 . 
In this paper, we describe a systematic procedure for de- 
scribing and investigating finite-area topological defects, 
and, through this method, identify several of these "new" 
defects as divacancies and divacancy complexes. 

The graphene structure can be represented as a hon- 
eycomb lattice (Fig. []Ja)). Every planar lattice has a 
topologically equivalent dim/[22], generated by convert- 
ing every n- vertex (vertex where n edges meet) to an 
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FIG. 1: (a) Ideal graphene honeycomb lattice with its dual 
tiling triangle tiling superimposed, (b) Graphene with Stone- 
Wales defect with its dual tiling superimposed, (c) Graphical 
representation of the Stone- Wales defect. 



n-tile (tile with n sides), and vice versa. Each edge in 
a tiling has exactly one corresponding edge in its dual. 
The dual of the honeycomb tiling is a regular triangu- 
lar tiling (Fig. [lja)). Similarly, every planar sp 2 bonded 
carbon structure (three bonds per carbon) has a dual lat- 
tice that is composed solely of triangles. The structure of 
topological defects in graphene can thus be represented 
in terms of how they change a regular triangular lattice 
into another triangular lattice. 

Consider for example, the Stone- Wales [23] defect 
(Fig. [TJb)). It rotates a graphene C-C bond by 90°. 
Its dual-lattice equivalent is the rotation of a patch of 
two joined triangles by 90° (Fig. [ljc)). The dual rela- 
tionship between two lattices is mutual. Therefore, it is 
possible to design a sp 2 -type defect structure in "dual 
space" by replacing a chosen patch of the ideal trian- 
gle tiling with a different triangulated patch of the same 
perimeter, and then taking its dual. Metaphorically, one 
cuts a patch out of graphene and then "grafts" or "trans- 
plants" a different sp 2 -bonded patch with the same num- 
ber of dangling bonds. Of particular interest is regions 
where the dual space replacement patch is also a portion 
of an ideal triangular lattice. Because the replacement 
structure is graphitic in this case, the formation energy 
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FIG. 2: (a) Patch of triangular tiling. (b;c) Two replacement 
patches, related by "most compatible donor" procedure, (d- 
f) Corresponding graphene structures. As described in the 
text, and seen here, the number of triangles surrounding a 
given vertex in dual space (the sum of the interior and exte- 
rior counts shown) equals the perimeter of the corresponding 
carbon ring in the defect structure. 



should be relatively low. In such cases, one effectively 
has a small grain of graphene inside bulk graphene with a 
closed grain boundary separating the regions. We called 
such defects "grain boundary loops" in Ref. \19\ 

There are clearly an infinite number of possible grain 
boundary loops in graphene. In dual space, they are de- 
scribed by the shapes of the (same-perimeter) triangular 
regions of the inner grain region "before" and "after" the 
procedure (Fig. [2]), with the relative orientation an addi- 
tional degree of freedom. Note that the "grafting" proce- 
dure provides a description of various topological defects 
and not how these defects actually form. This descrip- 
tion complements an alternative description of topologi- 
cal defects in terms of the sequence of elementary steps 
(divacancies, addimers, and Stone- Wales type bond ro- 
tations) that could generate a given defect out of ideal 
graphene [24]. 

Using the dual space formalism, it can easily be shown 
that the number of dangling bonds Ndb in a graphene 
patch satisfies Ndb = 3iV"c — 27V/5, where Nc is the num- 
ber of carbon atoms and Nf b is the number of "full" 
bonds in the interior. As a corollary to this formula, Nc 
must be even (odd) if Ndb is even (odd). Therefore, grain 
boundary loops, which preserve Ndb, can only change the 
number of carbon atoms by a multiple of 2. 

A n-vertex in dual space corresponds to a n-ring in 
graphene, that is, a ring of n bonded carbon atoms. 
The number of bonds emanating from a vertex in dual 
space equals the number of triangles surrounding that 
vertex. The number of triangles surrounding a vertex in 
the defect structure equals the number of exterior trian- 
gles sharing that vertex in the "before" state N e b plus 
the number of interior triangles sharing that vertex in 



the "after" state N ia (see Fig.|2|. Since N eb = 6 - N ib , 
where N ib is the number of "before" internal triangles, 
n = 6 - N ib + N ia . 

Studies of energies of various fullerenes containing 5- 
rings and 7-rings have led to the rules that the number 
of adjacent 5-5 pairs should be minimized and adjacent 
5-7 pairs maximized [25 . A similar empirical rule has 
been noted for topological defects in graphene [T9] [26] : 
it is energetically favorable for grain boundary loops to 
consist of alternating 5-rings and 7-rings. In terms of 
the above formula, alternating 5-rings and 7-rings are 
obtained if Ni a — Ni b is alternately +1 and -1 for the ver- 
tices on the perimeter of the replacement triangular patch 
in dual space. Continuing the transplant metaphor, we 
define "most compatible donor (s)" for each patch of the 
ideal triangle lattice (if any exist), as those patches whose 
perimeter vertices alternately span Ni a ± 1 and N{ a =p 1 
triangles with respect to the perimeter vertices of the 
original patch. Because either choice of sign is possible, a 
structure can have up to two inequivalent most compati- 
ble donors. Fig. [2] shows the two most compatible donors 
for one representative patch. The resultant defects are 
distinct; one reduces the number of carbon atoms by 2 
and the other reduces it by 4. 

We hypothesize that the most significant finite-area 
topological defects in graphene are grain boundary loops 
described by the "most compatible donor" procedure. 
We have enumerated all such defects with C2 V or higher 
symmetry that are contained within a diameter of 1.0 
nm or smaller, and that change the number of carbon 
atoms by two or fewer. We then investigated their stabil- 
ities and predicted their scanning tunneling microscopy 
(STM) signatures using density functional theory. 

First principles density functional theory (DFT) cal- 
culations, as encoded in the Vienna ab-initio Simu- 
lation Package (VASP) software [27]. were used per- 
form relaxations and to calculate the local density of 
states for STM simulations. Vanderbilt-type ultrasoft 
pseudopotentials[28] were used with a plane wave basis 
set with a cutoff energy of 211.1 eV. Calculations were 
mostly performed on supercells with approximately 216 
and 486 atoms. The results were extrapolated to esti- 
mate formation energies for isolated defects and their 
uncertainties. Efficient Brillouin zone integration was 
performed using a mesh containing six k-points in each 
supercell Brillouin zone ((1/3,0,0) and symmetry equiva- 
lents). These parameters were sufficient to obtain forma- 
tion energies in good agreement with previous results [29 . 
Further methodological details are given in [19]. 

The in-plane lattice parameter was set to the value for 
graphite, corresponding to epitaxial multilayer growth 
conditions. Structures with more atoms than ideal 
graphene were allowed to relax out-of-plane, while struc- 
tures with the same number of atoms or fewer atoms were 
fixed to remain planar. There are reports that such de- 
fects may prefer to relax out of plane [30 . However, weak 
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buckling instabilities might be suppressed under exper- 
imental conditions by adhesion of the surface layer in 
few-layer graphene to to the layer below. This is an issue 
for further exploration. 

STM images were simulated for comparison with ex- 
perimental fixed- voltage topographs. The tunneling cur- 
rent is approximated as proportional to the local elec- 
tronic density of states integrated between Ep and 
Ep + eAF, with Ep the Fermi level and AV the bias 
voltage[3T]. For these calculations, 4x4x1 grids of 
k-points in the Brillouin zone supercell were used, cen- 
tered at the origin. Using the experimental value of 
eAV = +0.3 eV, the best agreement with experimental 
STM images in [7 comes from setting Ep = Ep> +0.05 eV 
(Ep> the Dirac level) rather than the experimental value 
E F = Ep> + 0.3 eV. A significant part of the discrep- 
ancy comes from simulating a monolayer rather than the 
experimental case of Bernal-stacked few-layer graphene. 
Coupling between graphene layers in Bernal-stacked bi- 
layer graphene raises the energy of the low-lying anti- 
bonding states of the carbon atoms on the stacked sub- 
lattice by about 0.22 eV[32]. Aside from the offset dis- 
crepancy, very good agreement between the monolayer 
simulation and the few-layer experiment is found with 
greatly reduced computational cost compared with mul- 
tilayer DFT calculations. 

The various small high-symmetry graphene defects are 
shown in Fig. [3J The figures show the graphical rep- 
resentation of each defect, its relaxed atomic structure, 
the change in the number of atoms with respect to ideal 
graphene, and the formation energy. The structures that 
increase the number of atoms all have out-of-plane distor- 
tions of order 0.2 nm to 0.3 nm. Their relaxed structures 
have the same "rounded hill" appearance shown for such 
defects in the literature [14], and are not shown here. The 
corresponding STM simulations of the defects are shown 
in Fig. [4j The height scale is larger for those structures 
that distort out-of-plane. 

The defects shown side-by-side in Fig. [3] have graphical 
representations that are the inverse of each other, except 
for the Stone- Wales and flower defects, which are equiva- 
lent to their own inverses (as an aside, the addimer defect, 
frequently termed the "inverse Stone- Wales defect [H]", 
is not technically the inverse of the Stone- Wales defect, 
but the inverse of the divacancy ) . Using the dual formal- 
ism, it is easily seen that if a topological defect adds n 
atoms, its inverse removes n atoms. 

While the most compatible donor procedure (Fig. [2| 
was designed to yield structure with alternating 5 and 7 
rings, in the case of the divacancy, a patch of two triangles 
in dual space collapses into what is effectively a zero-area 
polygon with 4 sides (Fig. J3^c)). The resultant structure 
has a new polygon, namely an 8-ring, at the center. 

The most compatible donor procedure yields dual 
structure patches whose triangles are rotated 30° with 
respect to those of the original dual. When the patches 




FIG. 3: Small high-symmetry topological defects in graphene: 
(a) Stone- Wales [23]. (b) Flower [20]. (c) Divacancy. (d) Ad- 
dimer. (e) V 2 (555-777) [29] [33]. (f) H 5 , 6 ,r[34]. (g) V 2 (5555-6- 
7777) 29 . (h) 2-hexagon addimer reconstruction [35]. Defects 
related by inversion are side- by-side. Opposite and same-side 
triangle patch are the before-and-after representation in dual 
space. Next to the relaxed structures are the change in num- 
ber of atoms (above) and the formation energy (below), in eV, 
with uncertainty in units of 0.1 eV in parentheses (if greater 
than 0.05 eV). 

have sufficiently high symmetry, strain mismatch is min- 
imized when the transplant patch is oriented exactly 30° 
with respect to the surrounding structure. The procedure 
in this work thus explains the experimental observation 
of of 30° rotated regions inside grain boundary loops [26 . 

Two divacancy reconstructions, V<2(555 — 777) and 
^(5555 — 6 — 7777) are found that are lower energy than 
the simple divacancy. Such energy-lowering reconstruc- 
tions have previously been reported [29, 33 . In contrast 
to [29], the V<2(5555 — 6 — 7777) defect energy is lower 
than the ^(555 — 777) one for both periodic cells in- 
vestigated, although the difference between the extrapo- 
lated energies is less than their uncertainties. Since the 
energy of a grain boundary loop tends to increase with 
its perimeter [19], it is unlikely that there is any undis- 
covered divacancy reconstruction that has lower energy 
(unless it has lower symmetry than explored here). 

For an addition of two atoms, the addimer defect has 
lowest energy. Reconstructing the addimer gives config- 
urations that are higher in energy. Note, however, that 
both addimer reconstructions in Fig. [3] have appeared in 
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FIG. 4: (a-h) Simulated STM topographs for the defects 
shown in Fig. [3] (a-h), respectively. Height range from bright- 
est to darkest is 0.32 nm in nonplanar defects (d),(f), and (h), 
and 0.12 nm otherwise. All images in Fig. ^] to Fig. [7] show 
1.0 x 1.0 nm patches centered on a defect. 



numerical simulations of addimers on carbon nanotubes 
under strain [35]. 

Topological defects can be combined to make more 
complex defect structures. In the dual space descrip- 
tion nonoverlapping regions are simultaneously retrian- 
gulated. It is an open question whether and under 
which circumstances various topological defects attract 
or repel, but it is interesting that clustering has been 
observed experimentally for flower defects [18 and for 
divacancies[26]. 

Fig.[5]shows two defects observed in experimental STM 
topographs [7 that are well-matched by simulations of 
topological defects shown in Fig. [4] The sixfold defect 
shown in Fig. [5|a) agrees with the simulation of the 
flower defect (Fig. [5jb)), as previously shown in [19]. 




FIG. 5: (a) Experimental STM images of sixfold defect in 
graphene (all experimental images shown in Fig. [5}Fig. [7] are 
adapted from Ref. 7 ]). (b) Simulated STM image of flower 
defect is a good match for (a), (c) Experimental STM image 
of small twofold defect in graphene matches (d) simulated 
STM image of divacancy 
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FIG. 6: (a) Experimental STM image of large twofold de- 
fect in graphene. (b) The simulated STM image of a double 
divacancy agrees well with (a). 



The twofold symmetric defect in Fig. [5|c) agrees very 
well with the simulation of the divacancy (Fig. (5^d)). 
The characteristic dumbbell shape of the central region 
was previously predicted by Amara et al [36] , and distin- 
guishes the STM image of this defect from that of the 
twofold symmetric Stone- Wales defect. 

A second twofold defect observed in [7] is shown in 
Fig.^a). The double divacancy shown in Fig.[6jb) is an 
excellent match to the experimental image. The double 
divacancy has also been observed in single-layer graphene 
created by mechanical cleavage and then irradiated [26J. 
It contains a rectangle of 4 carbon atoms in its center. 
Note that the simulated STM image is not simply a su- 
perposition of the STM images of two individual diva- 
cancies. 

Finally, Fig. [7^a) presents the experimental STM im- 
age of a defect with threefold symmetry. This defect 
was initially modelled as a triple divacancy (Fig.[7^b)) in 
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FIG. 7: (a) Experimental STM image of threefold defect in 
graphene. (b) Simulated STM image of a triple divacancy 
matches central concave region well but matches poorly out- 
side this region, (c) Simulated local density of states of 
pseudo-threefold vacancy between Ef + 0.1 eV and Ef + 
0.2 eV matches peaks and nodes of experimental image very 
well 



analogy with the double divacancy. The triple divacancy 
has a triangle of carbon atoms in the middle. While a 
planar triangular carbon configuration is perhaps unex- 
pected on energetic grounds, a C-C-C triangle naturally 
occurs for the low-energy nonplanar bridging configura- 
tion of a single carbon adatom on graphene [37| [38] . and a 
C-C-C triangle also occurs in the cyclopropane molecule. 
Furthermore, the calculated formation energy of a triple 
divacancy in the configuration shown is 15.9 ± 0.5 eV, 
less than three times the formation energy of three iso- 
lated divacancies, 7.2 ± 0.2 eV each (Fig. |3|, giving a 
plausibility argument for the arrangement in Fig. (7^b). 
(The calculated energy of the double divacancy, 11.5 ± 
0.9 eV is also less than the total for the isolated divacan- 
cies, further showing the tendency of divacancies to clus- 
ter in graphene.) The simulated STM image in Fig. [7jb) 
reproduces the central concave triangle seen experimen- 
tally. The peaks and nodes outside this region, however, 
are poorly reproduced. 

An alternate explanation for the experimental three- 
fold defect is a single vacancy (Fig.^c)). Although a va- 
cancy in graphene relaxes in a way that breaks threefold 
symmetry [36, 39 , the calculated local density of states 
between + 0.1 eV and E^ + 0.2 eV is found to have 
pseudo-threefold symmetry (Fig. |7Jc)). (The DFT cal- 
culations for this simulation included spin polarization 
to incorporate magnetism.) The electronic structure in 
this energy range reproduces the concave central region 
and furthermore, the peak and node structure outside 
this region matches that of the experimental STM image. 
The simulated image in Fig. [T^c)) looks almost exactly 
like that simulated for a generic impurity on the surface 
of graphite by Mizes and Foster [40]. Their qualitative 
model was based on a tight-binding approximation and 
is equally valid for an adatom directly above a carbon, 
a vacancy, or a substitution. It is concluded that the 
threefold experimental defect is one of these three types 
that is located on a single carbon site. A problem with 
interpreting the defect as a vacancy is that, when sim- 
ulated to match the experimental bias of 0.3 eV, either 
from E£> + 0.05 eV and E^ + 0.35 eV as for the other 



STM simulations, or for any other 0.3 eV range near Ed, 
then the agreement with experiment is lost. A possible 
explanation for the discrepancy is strong finite-size unit 
cell artifacts, as the perturbation due to a vacancy falls 
off as l/rjH^ more slowly than for topological defects 
that preserve sp 2 binding. STM images for a single va- 
cancy in graphene need to be simulated for larger unit 
cells to test this hypothesis and compared with similar 
simulations for adatoms and substitutions. 

In conclusion, we introduce a powerful dual space 
method for describing finite-area topological defects in 
graphene. This method allows us to systematically ex- 
plore candidate low-energy defects and defect clusters. 
Previously unidentified defects in graphene were identi- 
fied as divacancies or divacancy clusters. 
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